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abstract 

Jets coexist with planetary scale waves in the turbulence of planetary atmospheres. The coherent compo¬ 
nent of these structures arises from cooperative interaction between the coherent structures and the incoherent 
small-scale turbulence in which they are embedded. It follows that theoretical understanding of the dynamics 
of jets and planetary scale waves requires adopting the perspective of statistical state dynamics (SSD) which 
comprises the dynamics of the interaction between coherent and incoherent components in the turbulent state. 
In this work the S3T implementation of SSD for baroti'opic beta-plane turbulence is used to develop a theory 
for the jet/wave coexistence regime by separating the coherent motions consisting of the zonal jets together 
with a selection of large-scale waves from the smaller scale motions which constitute the incoherent compo¬ 
nent. It is found that mean flow/turbulence interaction gives rise to jets that coexist with large-scale coherent 
waves in a synergistic manner. Large-scale waves that would exist only as damped modes in the laminar jet 
are found to be transformed into exponentially growing waves by interaction with the incoherent small scale 
turbulence which results in a change in the mode structure allowing the mode to tap the energy of the mean jet. 
This mechanism of destabilization differs fundamentally and serves to augment the more familiar S3T insta¬ 
bilities in which jets and waves arise from homogeneous turbulence with energy source exclusively from the 
incoherent eddy field and provides further insight into the cooperative dynamics of the jet/waves coexistence 
regime in planetary turbulence. 


1. Introduction 

A regime in which jets, planetary scale waves and vor¬ 
tices coexist is commonly observed in the turbulence of 
planetary atmospheres with the banded winds and embed¬ 
ded vortices of Jupiter and the Saturn north polar vortex 
constituting familiar examples (Vasavada and Showman 
2005; Sanchez-Lavega et al. 2014). Planetary scale waves 
in the jet stream and vortices such as cutoff lows are also 
commonly observed in the Earth’s atmosphere. Conserva¬ 
tion of energy and enstrophy in undamped 2D turbulence 
implies continual transfer of energy to the largest available 
spatial scales (Fjprtoft 1953). This upscale transfer pro¬ 
vides a conceptual basis for expecting the largest scales 
to become increasingly dominant as the energy of turbu¬ 
lence forced at smaller scale is continually transferred to 
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the larger scales. However, the observed large-scale stmc- 
ture in planetary atmospheres is dominated not by incoher¬ 
ent large-scale turbulent motion as would be expected to 
result from the incoherent phase relation of Fourier modes 
in a turbulent cascade, but rather by coherent zonal jets, 
vortices and waves of highly specific form. Moreover, 
the scale of these coherent structures is distinct from the 
largest scale permitted in the flow. An early attempt to un¬ 
derstand the formation of jets in planetary turbulence did 
not address the structure of the jet beyond attributing the 
jet scale to arrest of the incoherent upscale energy cascade 
at the length scale set by the value of the planetary vor- 
ticity gradient and a characteristic flow velocity (Rhines 
1975). In Rhines’s interpretation this is the scale at which 
the turbulent energy cascade is intercepted by the forma¬ 
tion of propagating Rossby waves. While this result pro¬ 
vides a conceptual basis for expecting zonal structures 
with spatial scale limited by the planetary vorticity gra¬ 
dient to form in beta-plane turbulence, the physical mech¬ 
anism of formation, the precise morphology of the coher- 
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ent structures and their stability are not addressed by these 
general considerations. 

Our goal in this work is to continue development of a 
general theory for the formation of finite amplitude struc¬ 
tures in planetary turbulence, specifically addressing the 
regime in which jets and planetary waves coexist. This 
theory identifies specific mechanisms responsible for for¬ 
mation and equilibration of coherent stmctures in plan¬ 
etary turbulence. A number of mechanisms have been 
previously advanced to account for jet, wave and vortex 
formation. One such mechanism which addresses exclu¬ 
sively jet formation is voiticity mixing by breaking Rossby 
waves leading to homogenization of potential vorticity 
(PV) in localized regions (Baldwin et al. 2007; Dritschel 
and McIntyre 2008) resulting in the case of barotropic 
beta-plane turbulence in broad retrograde parabolic jets 
and relatively narrow prograde jets with associated stair¬ 
case structure in the absolute vorticity. While PV stair¬ 
cases have been obtained in some numerical simulations 
of strong jets (Scott and Dritschel 2012), vorticity mixing 
in the case of weak to moderately strong jets is insufficient 
to produce a prominent staircase structure. Moreover, jets 
have been shown to form as a bifurcation from homoge¬ 
neous turbulence in which case the jet is perturbative in 
amplitude and wave breaking is not involved (Farrell and 
loannou 2003, 2007). 

Equilibrium statistical mechanics has also been ad¬ 
vanced to explain formation of coherent structures 
e.g. by Miller (1990) and Robert and Sommeria (1991). 
The principle is that dissipationless turbulence tends to 
produce configurations that maximize entropy while con¬ 
serving both energy and enstrophy. These maximum en¬ 
tropy configurations in beta-plane turbulence assume the 
coarse grained stmcture of zonal jets (cf. Bouchet and Ve- 
naille (2012)). However, the relevance of these results 
to the formation, equilibration and maintenance of jets in 
strongly forced and dissipated planetary flows remains to 
be established. 

Zonal jets and waves can also arise from modulational 
instability (Lorenz 1972; Gill 1974; Manfroi and Young 
1999; Berloff et al. 2009; Connaughton et al. 2010). This 
instability produces spectrally non-local transfer to the un¬ 
stable structure from forced waves and therefore presumes 
a continual source of waves with the required form. In 
baroclinic flows, baroclinic instability has been advanced 
as the source of these waves (Berloff et al. 2009). From the 
broader perspective of the statistical state dynamics the¬ 
ory used in this work, modulational instability is a special 
case of an S3T structural instability (Parker and Krommes 
2014; Parker 2014; Bakas et al. 2015). However, mod¬ 
ulational instability does not include the mechanisms for 
realistic equilibration of the instabilities at finite amplitude 
although a Landau-type term has been used to produce 
equilibration of the modulational instability (cf. Manfroi 
and Young (1999)). 


Another approach to understanding the jet/wave co¬ 
existence regime is based on the idea that jets and 
waves interact in a cooperative manner. Such a dynamic 
is suggested for example by observations of a promi¬ 
nent wavenumber-5 disturbance in the Southern Hemi¬ 
sphere (Salby 1982). Using azonally symmetric two-layer 
baroclinic model Cai and Mak (1990) demonstrated that 
storm track organization by a propagating planetary scale 
wave resulted in modulation in the distribution of synop¬ 
tic scale transients configured on average to maintain the 
organizing planetary scale wave. The symbiotic forcing 
by synoptic scale transients which on average maintains 
planetary scale waves was traced to barotropic interactions 
in the studies of Robinson (1991) and Qin and Robinson 
(1992). While diagnostics of simulations such as these are 
suggestive, comprehensive analysis of the essentially sta¬ 
tistical mechanism of the symbiotic regime requires ob¬ 
taining solutions of the statistical state dynamics underly¬ 
ing it and indeed the present work identifies an underly¬ 
ing statistical mechanism by which transients are system¬ 
atically organized by a planetary scale wave so as to on 
average support that planetary scale wave in a spectrally 
non-local manner. 

Stochastic structural stability theory (S3T) provides a 
statistical state dynamics (SSD) based theory accounting 
for the formation, equilibration and stability of coherent 
structures in turbulent flows. The underlying mechanism 
of jet and wave formation revealed by S3T is the spec¬ 
trally non-local interaction between the large-scale struc¬ 
ture and the small-scale turbulence (Farrell and loannou 
2003). S3T is a non-equilibrium statistical theory based 
on a closure comprising the nonlinear dynamics of the 
coherent large-scale structure together with the consis¬ 
tent second-order fluxes arising from the incoherent ed¬ 
dies. The S3T system is a cumulant expansion of the 
turbulent dynamics closed at second order (cf. Marston 
et al. (2008)), which has been shown to become asymp¬ 
totically exact for large-scale jet dynamics in turbulent 
flows in the limit of zero forcing and dissipation and in¬ 
finite separation between the time scales of evolution of 
the large-scale jets and the eddies (Bouchet et al. 2013; 
Tangarife 2015). S3T has been employed to understand 
the emergence and equilibration of zonal jets in planetary 
turbulence in barotropic flows on a beta-plane and on the 
sphere (Farrell and loannou 2003, 2007, 2009a; Marston 
et al. 2008; Srinivasan and Young 2012; Marston 2012; 
Constantinou et al. 2014; Bakas and loannou 2013b; To¬ 
bias and Marston 2013; Parker and Krommes 2014; Ait- 
Chaalal et al. 2016), in baroclinic two layer turbulence 
(Farrell and loannou 2008, 2009c) and in drift-wave turbu¬ 
lence in plasmas (Farrell and loannou 2009b; Parker and 
Krommes 2013). It has been used in order to study the 
emergence and equilibration of finite amplitude propagat¬ 
ing non-zonal structures in barotropic flows (Bakas and 
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loannou 2013a, 2014; Bakas et al. 2015) and the dynam¬ 
ics of blocking in two-layer baroclinic atmospheres (Bern¬ 
stein and Farrell 2010). It has also been used to study the 
role of coherent structures in the dynamics of the 3D tur¬ 
bulence of wall-bounded shear flows (Farrell and loannou 
2012; Thomas et al. 2014, 2015; Farrell et al. 2015). 

In certain cases a barotropic S3T homogeneous turbu¬ 
lent equilibrium undergoes a bifurcation in which non- 
zonal coherent structures emerge as a function of turbu¬ 
lence intensity prior to the emergence of zonal jets and 
when zonal jets emerge a new type of jet/wave equilib¬ 
rium forms (Bakas and loannou 2014). In this paper we 
use S3T to further examine the dynamics of the jet/wave 
coexistence regime in barotropic beta-plane turbulence. In 
order to probe the jet/wave/turbulence dynamics in more 
depth a separation is made between the coherent jets and 
large-scale waves and the smaller scale motions which are 
considered to constitute the incoherent turbulent compo¬ 
nent of the flow. This separation is accomplished using 
a dynamically consistent projection in Fourier space. By 
this means we show that jet states maintained by turbu¬ 
lence may be unstable to emergence of non-zonal traveling 
waves and trace these unstable eigenmodes to what would, 
in the absence of turbulent fluxes, have been damped wave 
modes of the mean jet. Thus we show that the cooperative 
dynamics between large-scale coherent and small-scale in¬ 
coherent motion is able to transform damped modes into 
unstable modes by altering the mode structure allowing it 
to tap the energy of the mean jet. 

In this work we also extend the S3T stability analysis of 
homogeneous equilibria (Farrell and loannou 2003, 2007; 
Srinivasan and Young 2012; Bakas and loannou 2013a,b, 
2014; Bakas et al. 2015) to the S3T stability of jet equilib¬ 
ria. We present new methods for the calculation of the S3T 
stability of jet equilibria, extending the work of Farrell and 
loannou (2003); Parker and Krommes (2014), which was 
limited to the study of the S3T stability of jets only with 
respect to zonal perturbations, to the S3T stability of jets 
to non-zonal perturbations. 

2. Formulation of S3T dynamics for barotropic 

j3 -plane turbulence 

Consider a non-divergent flow u = (m,v) on a /3-plane 
with coordinates x = {x,y)', in which x is the zonal direc¬ 
tion and y the meridional direction, and with the flow con¬ 
fined to a periodic channel of size InL x 2kL. The ve¬ 
locity field can be obtained from a streamfunction iff as 
u = z X Vv^, with z the unit vector normal to the {x,y) 
plane. The component of vorticity normal to the plane 
of motion is — dyii is given as = A\jf with 

A itef ^2 -i_ ^2 Laplacian operator. In the presence of 
dissipation and stochastic excitation the vorticity evolves 


according to: 

5,C = -u-VC-j3v-rC + vAC + v^^ , (1) 

in which the flow is damped by Rayleigh dissipation 
with coefficient r and viscous dissipation with coefficient 
V. The stochastic excitation maintaining the turbulence, 
^(x,r), is a Gaussian random process that is temporally 
delta-correlated with zero mean. 

Equation (1) is non-dimensionalized using length 
scale L and time scale T. The double periodic do¬ 
main becomes 2n x 2n and the non-dimensional vari¬ 
ables in (1) are: C* = u* = = 

e* = e/{L^T-'^), (5* = ^/{LT)-\ r* = 
rjT^^ and v* = V/{L?T^^), where asterisks denote non- 
dimensional units. Hereafter all variables are assumed 
non-dimensional and the asterisk is omitted. 

We review now the formulation of the S3T approxima¬ 
tion to the statistical state dynamics (SSD) of (1). The S3T 
dynamics was introduced in the matrix formulation by Far¬ 
rell and loannou (2003). Marston et al. (2008) showed that 
S3T comprises a canonical second-order closure of the ex¬ 
act statistical state dynamics and derived it alternatively 
using the Hopf formulation. Srinivasan and Young (2012) 
obtained a continuous formulation which facilitates ana¬ 
lytical explorations of S3T stability of turbulent statistical 
equilibria. 

An averaging operator by which mean quantities are ob¬ 
tained, denoted by angle brackets, (•), is required in order 
to form the S3T equations. Using this averaging operator 
the flow streamfunction is decomposed as 

C = z+C', (2) 

where 

Z(x,f) = (C(x,f)) , (3) 

is the mean field or the first cumulant of the vorticity and, 
similarly, for the derived flow fields, i.e. U, yf. The eddies, 

satisfy the important property that 

(C')=0, (4) 

which relies on the averaging operation satisfying the 
Reynolds condition (cf Ait-Chaalal et al. (2016)) that for 
any two fields / and g, 

{{f)8) = {f){8) ■ (5) 

The equation for the first cumulant is obtained by av¬ 
eraging equation (1), which after repeated use of (5) be¬ 
comes: 

d,Z+ IJ ■ VZ + pv + rZ- vAZ =-V • (u'C') , (6) 

in which we have we assumed (^} =0. The term 
— V • represents the source of mean vorticity aris¬ 

ing from the perturbation vorticity flux divergence. 
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The second cumulant of the vorticity fluctuations is the 
covariance 

C(x„,Xfo,f) = (C'(x„,f)C'(xi„f)) , (7) 

which is a function of five variables: time, t, and the coor¬ 
dinates of the two points x„ and xi,. We write (7) concisely 

as Cab = {CcXb)- 

All second moments of the velocities can he expressed 
as linear functions of C. For example the perturbation vor¬ 
ticity flux divergence source term, V • (u'(x,f)^'(x,f)), in 
the mean vorticity equation (6) can be written as a function 
of C as follows: 

V.{u'i;')=^-V.{u'aa+uiOa=b 

= ^-v.[zx{VaA-^+VbA,^){c:.a)]a.b 

= )^‘‘b]a=b 

=‘^ R{C) ( 8 ) 

def 

in which u'- = u'(xj,r) and subscripts in the differential 
operators indicate the specific independent spatial vari¬ 
able the operator is defined on. To derive (8) we used 
that u'- = z X ^jAj^^j, with the inverse Laplacian. 
The notation a = b indicates that the function of the five 
independent variables, x^, Xb and t, in (8) is to be consid¬ 
ered a function of two independent spatial variables and t 
by setting Xg =Xb = x. By denoting the divergence of the 
mean of the perturbation vorticity flux V • in (8) as 
R{C) we underline that the forcing of the mean vorticity 
equation (6) by the perturbations depends on the second 
cumulant (the covariance of the vorticity field). Adopting 
this notation for the divergence of the mean of the pertur¬ 
bation vorticity flux, the equation for the mean vorticity 
(the first cumulant) (6) takes the form: 

drZ + \J-VZ + pV + rZ-vAZ = -R{C) . (9a) 

The equation for the perturbation vorticity is obtained by 
subtracting (6) from (1): 

drC’ = - (U • VC' +u' . VZ) - V . (u'C' - (u'C')) 

- j3v' - rC' + vAC' -I- 

= AC'-V.(u'C'-(u'C'))+v'^^ , (9b) 

where 

A =-U-V- [/35^-(AU)-V]A-'-r-hvA. (10) 

Using (9b), definition (7) and noting that (C') =0 we ob¬ 
tain the evolution equation for C: 

drC={Cdra + adrQ 

= (Aq +Ab)C + Vs {Ca^b + Cb^a) 

+{[^a-{<Q]Cb+[^b-Wb)]Q ■ 

(11) 


Both terms {[V „ • (KCa)] Q and {[Vb-{KVb)]Q 
in (11) can be expressed as linear functions of the third 

cumulant of the vorticity fluctuations, Tgbc '= (ciac'), 

e-g- 

([v„.(u'c')]a) = 

= ^(v„-[u'c+u'c'],^,a) 

= \Va-[ix{VaAV+VaAV){CaO]c^a 

= i V„ . [z X {VaA- ' + V,AC' ) Fabc] a^a > ( ^2) 

explicitly revealing that the dynamics of the second cu¬ 
mulant of the vorticity fluctuations is not closed. Nota¬ 
tion c ^ a indicates that the function of independent spa¬ 
tial variables Xg, Xb and x^ should be considered a func¬ 
tion of only Xa and Xb after setting Xc —)• x^. The S3T 
system is obtained by truncating the cumulant expansion 
at second order either by setting the third cumulant term 
in (11) equal to zero or by assuming that the third cumu¬ 
lant term is proportional to a state independent covariance 
Q{xa,Xb). The latter is equivalent to representing both the 
nonlinearity, V • (u'C' — (u'C'))> externally im¬ 

posed stochastic excitation together as a single stochastic 
excitation V^^{x,t) with zero mean and two point and 
two time correlation function: 

(C(x„,fl)C(xi„f2)) = S{ti-t2)Q{Xa,Xb) , (13) 

from which it can be shown that^: 

(C^b + Vb^a) = ViQVa.^b) , ( 14 ) 

and consequently (11) simplifies to the time dependent 
Lyapunov equation: 

5rC=(A,+Afc)C + ee, (15) 

in which the subscripts a, b on C and Q are implied. 

Using parametrization (13) to account for both the per¬ 
turbation nonlinearity, V • (u'C' — (u'C')), and the exter¬ 
nal stochastic excitation, implies that full correspon¬ 
dence between the mean equation (9a) coupled with the 
parameterized perturbation equation (9b) and the nonlin¬ 
ear dynamics (1) requires that the stochastic term accounts 
fully for modification of the perturbation spectrum by the 
perturbation nonlinearity in addition to the explicit exter¬ 
nally imposed stochastic excitation. It follows that the 
stochastic parameterization required to obtain agreement 
between the approximate statistical state dynamics and 
the nonlinear simulations differs from the explicit external 
forcing alone unless the eddy-eddy interactions are negli¬ 
gible. 


'Assumption (13) implies identity (14) even when C obeys the non¬ 
linear (9b) (cf. Farrell and loannou (2014); Constantinou (2015). 
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The resulting S3T system is an autonomous dynamical 
system involving only the first two cumulants that deter¬ 
mines their consistent evolution. The S3T system for a 
chosen averaging operator is; 

dtZ=-\J-VZ-l5V-rZ+vAZ-R{C) , (16a) 

d,C={A,+Ai,)C + eQ. (16b) 

For the purpose of studying turbulence dynamics it is 
appropriate to choose an averaging operator that isolates 
the physical mechanism of interest. Typically the averag¬ 
ing operator is chosen to separate the coherent structures 
from the incoherent turbulent motions. Coherent struc¬ 
tures are critical components of turbulence in shear flow 
both in the energetics of interaction between the large and 
small scales and in the mechanism by which the statisti¬ 
cal steady state is determined. Retaining the nonlinear¬ 
ity and structure of these flow components is crucial to 
constructing a theory of shear flow turbulence that prop¬ 
erly accounts for the role of the coherent structures. In 
contrast, nonlinearity and detailed structure information is 
not required to account for the role of the incoherent mo¬ 
tions and the statistical information contained in the sec¬ 
ond cumulant suffices to include the influence of these on 
the turbulence dynamics. This results in a great practical 
as well as conceptual simplification that allows a theory 
of turbulence to be constructed. In the case of beta-plane 
turbulence a phenomenon of interest is the formation of 
coherent zonal jets from the background of incoherent tur¬ 
bulence. To isolate the dynamics of jet formation zonal av¬ 
eraging is appropriate. Alternatively, if the focus of study 
is the emergence of large planetary scale waves the av¬ 
eraging operation would be an appropriate extension of 
the Reynolds average over an intermediate spatial scale to 
produce a spatially coarse grained /fine grained flow sep¬ 
aration. An averaging operation of this form was used by 
(Bernstein and Farrell 2010) in their S3T study of block¬ 
ing in a two-layer baroclinic atmosphere and by Bakas 
and loannou (2013a, 2014) to provide an explanation for 
the emergence of travelling wave structures (“zonons”) in 
barotropic turbulence. However, the Reynolds average de¬ 
fined over an intermediate time or spatial scale; 

21 Jt-T 
or 

Hof 1 rx+X ry+y 

(/(X,f)) = ^ ^ Jy-Y 

satisfies the Reynolds condition (5) only approximately 
and to the extent that there is adequate scale separation. 
The S3T system that was derived in (16) is exact if the av¬ 
eraging operation is the zonal average and an adequate ap¬ 
proximation for jets and a selection of large-scale waves if 


there is sufficient scale separation to satisfy the Reynolds 
condition (5). 

Because the scale separation assumed in (16) is only 
approximately satisfied in many cases of interest, an alter¬ 
native formulation of S3T will now be obtained in which 
separation into two independent interacting components 
of different scales is implemented (a similar formulation 
is used by Marston et al. (2016)). This formulation makes 
more precise the dynamics of the coherent jet and wave 
interacting with incoherent turbulence regime in S3T. 

The required separation is obtained by projecting the 
dynamics (1) on two distinct sets of Fourier harmonics. 
Consider the Fourier expansion of the streamfunction, 

= , ( 18 ) 

kx ky 

with k = {kx, ky) and the projection operator Pk defined as 
(cf. Frisch (1995)); 

PkV= Z (19) 

\kx\<K ky 

SO that the large-scale flow is identified through stream- 
function 'P = and the small-scale flow through \f/' = 
{I — Pk)^/ where; 

{I-Pk)W= Z (20) 

ky 

with I the identity. Similarly, vorticity and velocity fields 
are decomposed into =Z+l^' and u = U-I-u'. 

From (1) and under the assumption that the stochas¬ 
tic excitation projects only on the small scales, the large 
scales evolve according to; 

5,Z = -P;,[U.VZ + V.(u'C')] 

-P;f(U-VC' + u'-VZ) 
-PV-rZ+vAZ , (21a) 

while the small scales evolve according to; 

d,C'=-iI-PK){V-V^' + u'-VZ) 

-{I-Pk)[V-VZ+V-{u'C')] 
-Pv'-rC' + vAC' + V^^ . (21b) 

If Pk were an averaging operator that satisfied the 
Reynolds condition (5) term P/f (U • -fu'• VZ) 
in (21a) would vanish. Here it does not, as both of 
these terms scatter energy to the large scales. How¬ 
ever, an energetically closed S3T system for the first 
two cumulants can be derived by making the quasi- 
linear approximation (QL) in (21b), i.e. neglect the terms 
(/ — Pjf) [U • VZ-I- V • (u'l^')] that represent projection of 
the eddy-eddy and large-scale-large-scale interactions on 
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the eddy flow components, and additionally neglect the 
terms Pk{V • + u' • VZ) in the large-scale equa¬ 

tion (21a). These later terms as well as (I — Pk) (U • VZ) 
are not of primary importance to the dynamics and in any 
case vanish with sufficient scale separation. With these 
terms neglected we obtain the projected QL system: 

drZ = -pV-rZ+vAZ-PK[V-VZ+V- (u'C')] , 

(22a) 

drC'=-{l-PK) (U-VC' + u'-VZ) 

- j3v' - rC' -f vAC' -f , (22b) 

which conserves total energy and enstrophy in the absence 
of forcing and dissipation. The conservation properties of 
the full barotropic equations are retained because the typ¬ 
ically small terms that have been discarded scatter energy 
and enstrophy between (21). 

Assuming Z = Pk{Q is the coherent flow and C = 
the covariance of the incoherent eddies, 
with ( • ) being an average over forcing realizations, we 
obtain the corresponding S3T system for the first two cu- 
mulants: 

d,Z = -pV-rZ + vAZ-PK[lJ-VZ + R{C)] , (23a) 

d,C = (/ - PKa)Aa C + (/ - PKb)Ab C + eQ. (23b) 

It can be shown that these equations have the same 
quadratic conservation properties as the S3T equa¬ 
tions (16) and the full nonlinear equations (1). Note 
that for = 0 this projection formulation reduces to the 
zonal mean/perturbation formulation employed previously 
to study zonal jet formation (Farrell and loannou 2003, 
2007; Srinivasan and Young 2012). 


3. Specification of the parameters used in this work 

Assume that the large-scale phase coherent motions 
occupy zonal wavenumbers \kx\ = 0,1 and all zonal 
wavenumbers \kjc\ > 2 represent phase incoherent motions, 
so that Pk has K = 1. 

The covariance of the stochastic excitation 
in (13) is assumed to be spatially homogeneous, 
i.e. Q{xa,Xb) = Q(xa —Xb), and can be associated with its 
Fourier power spectrum 0(k): 

G(x, - x&) = I ^ . (24) 


Unless otherwise indicated calculations are performed 
with the anisotropic power spectrum: 


G(k) 


{47t/Nf)k^e-'‘^‘‘^ 

-erf(M) 


L 

kf^Kf 


\5{kx — kf) + 5{kx + kf)^ , 
(25) 


with k = |k|, d = 0.2, Kf = {2,3,..., 14} the zonal 
wavenumbers that are forced and Nf the total number of 


excited zonal wavenumbers. This spectrum is biased to¬ 
wards small ky numbers consistent with the assumption 
that the forcing arises from baroclinic growth processes. 
The spatial excitation covariance, Q(k), has been normal¬ 
ized so that each kf injects equal energy and the total en¬ 
ergy injection rate is unity, i.e. 2(k) satisfies^: 


f 6(k) 

J (27r)2 2k^ 


(26) 


With this normalization the rate of energy injection by the 
stochastic forcing in (1), (16), (21), (22) and (23) is e and 
is independent of the state of the system because ^ has 
been assumed temporally delta-correlated. 

We choose j3 = 10, r = 0.15 and v = 0.01 as our pa¬ 
rameters. For L = 1200km and T = 6d these correspond 
to /3 = 1.6 X 10^^’ m^'s^^ and an e-folding time for lin¬ 
ear damping of 40d. The diffusion coefficient v = 0.01 is 
chosen so that scales of the order of the grid are damped 
in one non-dimensional time and it corresponds to an e- 
folding time for scales of the order of 1000 km (non- 
dimensional wavenumber k^ = 1 in our channel) of ap¬ 
proximately 400 d. With these parameters the channel has 
zonal extent about 7500 km, which corresponds to 1/4 of 
the latitude circle at 45°, one unit of velocity corresponds 
to 23 ms^* and non-dimensional e = 1 corresponds to an 
energy input rate of 1.03 x 10^^ Wkg^/ Simulations pre¬ 
sented in this work are performed using a pseudospectral 
code with = Ny = 64 grid points. 


4. S3T jet equilibria 

Fixed points of the S3T system correspond to statistical 
equilibria of the barotropic dynamics. We study these sta¬ 
tistical equilibria as a function of e. For all values of e and 
all homogeneous stochastic forcings there exist equilibria 
that are homogeneous (both in x and y) with 

fjft-fQQ-) rhr _x,)-p f ^ __ Jk-{Xa-Xi,) 

V -(U,U),c (X, Xb)-ej (2K)^2{r + vk^f 

(27) 

where Q(k) is the power spectrum of the stochastic forc¬ 
ing, defined in (24). 

However, these equilibria become unstable when e ex¬ 
ceeds a critical value. For values of e exceeding this 
critical value zonal jets arise from a supercritical bifur¬ 
cation (Farrell and loannou 2003, 2007; Srinivasan and 
Young 2012; Parker and Krommes 2013, 2014; Constanti- 
nou et al. 2014). These jets are constrained by the pe¬ 
riodic domain of our simulations to take discrete values 


stochastic term \/e^ with spatial covariance given by (13) can 
be shown to inject average energy per unit area in the fluid at a 
rate (L,-Ly)“*/d^x (l/rvYi') = e[(2;r)“^J d^k Q{k)/(2U)]. Since 
dimensional ^ has units we obtain from (13) that Q has 

dimensions therefore its Fourier transform Q is dimensionless. 
Hence (26) is valid for all values of the dimensional parameters. 
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Fig. 1. (a) Normalized turbulent energy input rates, e/fic, at which 
the homogeneous state becomes unstable to jet (tiji = 0) perturbations 
as a function of the jet meridional wavenumber riy. Dots indicate 
wavenumbers allowed in the model, Ec is the minimum energy input rate 
for jet emergence. Jets first emerge in an unrestricted eigencalculation 
at Ec = 0.2075 with unallowed wavenumber tiy = 2.82. For e/e^ <1.18 
the homogeneous state is stable to = 2 mean flow perturbations and 
iiy = 2 jet equilibria do not exist, (b) The riy = 2 zonal jet S3T equi¬ 
librium structure at e/Ec = 1.2,2,5,9,13.65 (marked with x in panel 
(a)). Increasing supercriticality results in increasing equilibrium jet am¬ 
plitude and deviation of the jet structure from the sinusoidal eigenmode 
form. 


of meridional wavenumber, riy. The critical curve in the 
{£,ny) plane separating the region in which only stable 
homogeneous turbulence equilibria exist from the region 
in which stable or unstable jet equilibria exist is shown 
for the chosen parameters in Fig. 1. This marginal curve 
was calculated using the eigenvalue relation for inhomo¬ 
geneous perturbations to the homogeneous S3T equilib¬ 
rium in the presence of diffusive dissipation, in the man¬ 
ner of Srinivasan and Young (2012) and Bakas and loan- 
nou (2014) with the wavenumber taking continuous 
values, but with the understanding that only integer val¬ 
ues of tiy satisfy the quantization conditions of the chan¬ 
nel. S3T instability of the homogeneous state first oc¬ 
curs at Hy = 2.82 for = 0.2075, which corresponds to 
2.15 X lO^^Wkg^*. Jets with rty = 3 emerge at 1.005ec 
and jets with tiy = 2 at I.ISEc. Examples of riy = 2 jet 
equilibria are shown in Fig. lb. The riy = 2 jet equilib¬ 
ria have mean flows and covariances that are periodic in y 
with period a = 7C and satisfy the time-independent S3T 


equations: 

I ') r] = rW - vd^W , (28a) 

{Al+Al)C = -eQ, (28b) 

with 

A" = -Wd, - [j3 - {d^U^)] - r + vA . (29) 

A basic property of the jet equilibria, which is shared by 
all S3T equilibria, is that they are hydrodynamically sta¬ 
ble (cf. Farrell and loannou (2014)). Stability is enforced 
at the discrete wavenumbers consistent with the finite do¬ 
main of the problem and not necessarily on the continuum 
of wavenumbers appropriate for an unbounded domain. 

5. S3T stability of the jet equilibria 

We are interested in the S3T stability of these = 2 
jet equilibria to non-zonal perturbations. The stability 
of jet equilibria to homogeneous in x perturbations has 
been investigated previously by Farrell and loannou (2003, 
2007) for periodic domains and by Parker and Krommes 
(2014); Parker (2014) for infinite domains. A comprehen¬ 
sive methodology for determining the stability of jet equi¬ 
libria to zonal and non-zonal perturbations was developed 
by Constantinou (2015). Recalling these results, perturba¬ 
tions (5Z, 5C) about the equilibrium state satis¬ 

fying equations (28), evolve according to^: 

d,8Z = PK[A'^ 8Z + R{8C)\ , (30a) 

d, 8C={I- Pko) {K8C + 8A,C) 

+ iI-PKb)iAl8C+8AbC) , (30b) 

with R as in (8)), A® defined in (29) and 

5A =^-5U-V-b(A5U)-VA-', (31) 

where 5U = z x VA^'5Z is the perturbation velocity 
field. 

Because of the homogeneity of the jet equilibria in the 
zonal, X, direction the mean flow eigenfunctions are har¬ 
monic functions in x and also because the equilibrium 
mean flow and covariance are periodic in y with period a, 
i.e. U‘^{y + a} = U‘^{y), Bloch’s theorem requires that each 
eigenfunction is a plane wave in y, modulated by a 
periodic function with period a in y (Cross and Greenside 
2009). Therefore, the eigenfunctions take the form: 

= (32a) 

— ^injc{xa-\-xi,)/2^iqyiya-\-yb)/2^(yt 

~ -^biyaiyt) ~ ^aiyb'^ya)'] 7 

(32b) 


^These perturbations equations are valid for equilibria inhomoge¬ 
neous in both X and v directions. In the case of jet equilibria the projec¬ 
tion operators are redundant. 
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Fig. 2. (a) Maximum S3T growth rates a,- as a function of e/Cc for ri.x — 0 and rix — I perturbations to the Uy — 2 equilibrium jets. The jet 
is unstable to jet {Ux — 0) perturbations for 1.18 < e/Ec ^ i-44 and e/Ec > 10.14 and to Ux — 1 wave perturbations for e/Ec > 6.80. (b) The 
corresponding phase speeds Cr of the m 
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Fig. 3. (a): Contour plot of the streamfunction of the least stable non-zonal rix — 1 S3T mean flow wave eigenfunction of the riy — 2 jet 
equilibrium at e/Ec — 1.2. This wave has growth rate G,- — —0.047 and phase speed Cr — —0.98. The equilibrium jet is shown in solid white. 
Positive (negative) contours are indicated with solid (dashed) lines and the zero contour with a thick solid line, (b): The power spectrum of 
the mean flow eigenfunction. This jet equilibrium is unstable to = 0 perturbations but stable to = 1 perturbations. The least stable rix — 1 
eigenfunction is Bloch qy — I with power at riy — 3. (c,d): Same for the least stable fix — 1 S3T eigenfunction of the equilibrium atE/Ec — 5. The jet 
is stable both totix — O and rix — 1 perturbations and the least stable rix — 1 eigenfunction (a,- = —0.033, Cr — —2.18) is Bloch qy — Q with power at 
fly — 2. (e,f): Same for the maximally growing rix — 1 S3T eigenfunction of the jet dXE/Ec — 9. The jet is stable to = 0 perturbations but unstable 
to = 1 perturbations and the most unstable fix — 1 eigenfunction (<7^ = 0.099, Cr — —3.81) is Bloch qy — I with power dXriy — 1. (g,h): Same for 
the maximally growing fix — ^ S3T eigenfunction of a strong equilibrium jet dXE/Ec — 13.65. The most unstable fix — i eigenfunction {Gr — 0.083, 
Cr — —5.99) is Bloch qy — 1 with power at fiy — 1. In this case, fix — 0 perturbations are more unstable with Gr — 0.324. 
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with \nx\ < K, 5Zn^^g^{y) periodic in y with period a and 
5Cn,,gy{xa -Xh,ya,yb) periodic in^a andy^, with period a. 
We have chosen 5C to he a symmetric function of x un¬ 
der the exchange Xa o The zonal wavenumber, n^, 
takes integer values in order to satisfy the periodic bound¬ 
ary conditions in x, and the Bloch wavenumber, qy, takes 
integer values in the interval < 7t/a, in order to sat¬ 
isfy the periodic boundary conditions in y (Constantinou 
2015). The eigenvalue a determines the S3T stability of 
the jet as a function of and qy. The jet is unstable when 

Or = Re(a) > 0 and the S3T eigenfunction propagates in 

def 

X with phase velocity = — Im(a)/«j: for «;c 7 ^ 0. When 
ny = 0 the eigenfunctions are homogeneous in the zonal, x, 
direction and correspond to a perturbation zonal jet. When 
tiyy^O the eigenfunctions are inhomogeneous in both and 
y and correspond to a wave. These perturbations, when un¬ 
stable, can form non-zonal large-scale structures that co¬ 
exist with the mean flow, as in Bakas and loannou (2014). 
For jets with meridional periodicity a = K, qy can take 
only the values = 0,1 and because these jets have a 
Fourier spectrum with power only at the even wavenum¬ 
bers, a qy = 0 Bloch eigenfunction has power only at even 
wavenumbers, while a = 1 Bloch eigenfunction has 
power only at odd wavenumbers. 

The maximum growth rate, Or, of the S3T eigenfunc¬ 
tion perturbations to the S3T equilibrium jet with tiy = 2 
(cf. Fig. lb) is plotted in Fig. 2a as a function of super¬ 
criticality e/Cc for both perturbations of jet form {tiy = 0) 
and non-zonal form (with n,^= 1). Consider first the sta¬ 
bility of the S3T jet to jet perturbation, that is to = 0 
eigenfunctions. Recall that the jets with Uy = 2 emerge at 
e/Cc = 1.18, and for e/Cc < 1.18 (shaded region in Fig. 2a) 
there are no ny = 2 equilibria. The dashed line shows 
the smallest decay/fastest growth rate of perturbations to 
the homogeneous equilibrium state that exists prior to jet 
formation. The most unstable eigenfunctions of the ho¬ 
mogeneous equilibria at these e are jets with wavenumber 
«y = 3 (not shown, cf. Fig. la). The small amplitude equi¬ 
librated «y = 2 jets that form when e marginally exceeds 
the critical e/Cc = 1.18 are unstable to jet formation at 
wavenumber = 3 with jet eigenfunction similar to the 
maximally growing «y = 3 eigenfunction of the homoge¬ 
neous equilibrium. This S3T instability of the small am¬ 
plitude tiy = 2 jet equilibria to «y = 3 jet eigenfunctions, 
which is induced by the «y = 3 instability of the nearby 
homogeneous equilibrium, was identified by Parker and 
Krommes (2014) as the universal Eckhaus instability of 
the equilibria that form near a supercritical bifurcation. 
The Eckhaus unstable S3T «y = 2 jets are attracted to the 
S3T «y = 3 stable jet equilibrium over the small interval 


^The covariance eigenfunction does not need to be symmetric or 
Hermitian in its matrix representation, but both symmetric and asym¬ 
metric parts have the same growth rate. For a discussion of the proper¬ 
ties of covariance eigenvalue problems cf. Farrell and loannou (2002). 


1.18 < e/Cc < 1.44. At higher supercriticalities in the in¬ 
terval 1.44 < e/Sc < 10.14 the tiy = 2 jets become sta¬ 
ble^ to tlx = 0 eigenfunctions. The jets eventually become 
unstable to = 0 eigenfunctions for e/Sc > 10.14. The 
most unstable «jc = 0 eigenfunction at e/Cc = 11 is a Bloch 
^y = 1 eigenfunction, dominated by a «y = 1 jet that will 
make the jets of the tiy = 2 equilibrium merge to form a 
«y = 1 jet equilibrium (cf. Earrell and loannou (2007)). 

The maximum growth rate of the jet equilibria to = 1 
non-zonal eigenfunctions is also shown in Eig. 2a. Unlike 
the jet eigenfunctions, which are stationary with respect to 
the mean flow, these eigenfunctions propagate retrograde 
with respect to the jet minimum; the phase velocity of the 
eigenfunction with maximum real part eigenvalue is plot¬ 
ted as a function of e/Cc in Eig. 2b. Eigenfunctions with 
tlx = 1 are stable for jets with e/Sc < 6.80 and when they 
become unstable the jet is still stable to jet {iix = 0) pertur¬ 
bations. The structure of the least damped/fastest grow¬ 
ing eigenfunctions at various e/Sc are shown in Eig. 3. 
In Eig. 3a,b is shown the least stable eigenfunction of the 
weak jet at e/Cc = 1.2. The eigenfunction is Bloch qy = l 
with almost all power at tiy = 3. The phase velocity of this 
eigenfunction is c,- = —0.98 which is slightly slower than 
the Rossby phase speed —1 (i.e., with j3 = 10, 

kx = l,ky = 3). In Eig. 3c,d is shown the least stable tix = l 
mode for the jet at e/Sc = 5 which is Bloch qy = 0 with 
almost all power at tiy = 2 and phase speed Cr = —2.18, 
which corresponds to a slightly modified Rossby phase 
speed with effective PV gradient of jSgff =10.9 instead of 
the j3 = 10 of the uniform flow. In Eig. 3e-h are shown the 
maximally unstable tix = I eigenfunctions for the jets at 
elec=9 and e/Ec = 13.65. Both are Bloch qy = 1 with al¬ 
most all power at tiy = 1. At e/Ec = 9 the mode is trapped 
in the retrograde jet, a region of reduced PV gradient, and 
the structure of this mode as well as its phase speed cor¬ 
responds, as shown in the next section, to that of an ex¬ 
ternal Rossby wave confined in this equilibrium flow. At 
e/Ec = 13.65 the eigenfunction is trapped in the prograde 
jet, a region of high PV gradient, and the stmcture of this 
mode as well as its phase speed corresponds to that of an 
external Rossby wave in this equilibrium flow. 

6. The mechanism destahilizing S3T jets to = 1 non- 

zonal perturbations 

We now examine the stability properties of the tiy = 2 
equilibrium jet maintained in S3T at e/Ec = 9. At e/Ec = 9 
the jet is stable to = 0 jet S3T perturbations but un¬ 
stable to tlx = I non-zonal perturbations with maximally 


^The periodic boundary conditions always allow the existence of 
a jet eigenfunction with zero growth and with structure that of the y 
derivative of the equilibrium jet and covariance. This eigenfunction 
leads to a translation of the equilibrium jet and its associated covari¬ 
ance in the y direction. This existence of this neutral eigenfunction can 
be a verified by taking the y derivative of (28). We do not include this 
obvious neutral eigenfunction in the stability analysis. 
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growing eigenfunction growth rate Or = 0.099 and phase 
speed Cr = —3.806, which is retrograde at speed 1.61 with 
respect to the minimum velocity of the jet. 

Because the jet U‘^ is an S3T equilibrium the operator A'' 
is necessarily stable to perturbation at zonal wavenumbers 
that are retained in the perturbation dynamics, i.e. e Kf. 
The maximum growth rate of operator A*^ as a function of 
kx for the jet e/Sc = 9 is shown in Fig. 4a, with the inte¬ 
ger valued wavenumbers that are included in the S3T dy¬ 
namics and are responsible for the stabilization of the jet 
indicated with a circle in this figure. This equilibrium jet, 
despite its robust hydrodynamic stability at all wavenum¬ 
bers, in both the mean and perturbation equations, and es¬ 
pecially its hydrodynamic stability to = 1 perturbations, 
is nevertheless S3T unstable at = 1. 

Although it is not formed as a result of a traditional hy¬ 
drodynamic instability, this S3T instability is very close in 
structure to the least stable eigenfunction of A'^ at = 1, 
as it can be seen in Fig. 4c,d. The spectrum of A‘^ at «^ = 1 
is shown in Fig. 4b. The eigenfunctions associated with 
this spectrum consist of viscous shear modes with phase 
speeds within the flow and a discrete number of external 
Rossby waves with phase speeds retrograde with respect 
to the minimum velocity of the flow (cf. Kasahara (1980)). 
In this case there are exactly 5 external Rossby waves with 
phase speeds Cr = —3.70,—9.80,—5.92,—2.33,—2.37 all 
decaying with k^Ci = —0.15,—0.16,—0.17,—0.18,—0.24 
respectively. We identify the S3T = 1 unstable eigen¬ 
function, shown in Fig. 4d, which has phase speed = 
—3.81 with S3T destabilization of the least stable of the 
external Rossby waves, shown in Fig. 4c, that has phase 
speed Cr = —3.70. This instability arises by Reynolds 
stress feedback that exploits the least damped mode of 
A‘^, which is already extracting some energy from the 
jet through the hydrodynamic instability process, thereby 
making it S3T unstable. This feedback process transforms 
a mode of the system that while extracting energy from 
the mean nevertheless was decaying at a rate k^ct = —0.15 
into an unstable mode growing at rate Or = 0.099. Consis¬ 
tently, note in Fig. 4d that the streamfunction of the S3T 
eigenfunction is tilting against the shear indicative of its 
gaining energy from the mean flow. 

We quantify the energetics of the S3T instability in or¬ 
der to examine the instability mechanism in more detail. 
The contribution to the growth rate of this 11^=1 eigen¬ 
function from interaction with the mean equilibrium jet is 

1 (Ai„v(t/«)5Z,5Z) + (5Z,Ai„v(C/«)5Z) 

^^“=2 - 

where (/, g) =^ —(27r)^^/d^x is the inner prod¬ 

uct in energy metric and 

Ai„v(C/) = -Udr, -[P- (d^U)] a,A-‘ , (34) 


is the inviscid part of (10) with V = 0. The contribution to 
the growth rate of the = 1 eigenfunction from Reynolds 
stress mediated interaction with the small scales is 


1 {5Z,R{dC)) + {R{dC),5Z) 
~2 {5Z,5Z) 


(35) 


The net growth rate of the perturbation = I eigenfunc¬ 
tion is then Or = CJio + CJi> -I- CTid, with 


1 (Ad 5Z, dZ) -k (5 Z,Ad5Z) 
“ 2 (dZ,dZ) 


(36) 


the loss to dissipation, where 


AD = -r+vA, (37) 

is the dissipation part of operator (10). 

For the S3T unstable eigenfunction shown in Fig. 4d, 
the growth rate Or = 0.099 arises solely from interac¬ 
tion with the mean flow, which contributes aio = 0.303, 
while the energy transfer from the small scale perturba¬ 
tion field contributes negatively, ai> = —0.016, with dis¬ 
sipation accounting for the remainder aio = —0.188. In¬ 
terestingly, this S3T unstable mode is solely supported in 
its energetics by induced non-normal interaction with the 
mean jet and loses energy to the Reynolds stress feedback 
which is responsible for the instability. This remarkable 
mechanism arises from eddy flux interaction transform¬ 
ing damped waves into exponentially growing waves by 
changing the wave structure so as to tap the energy of 
the mean jet. This novel mechanism destabilizes the wave 
even though the direct effect of the Reynolds stresses is to 
stabilize it. This mechanism of destabilization differs from 
that acting in more familiar S3T instabilities in which jets 
and waves arise directly from their interaction with the in¬ 
coherent eddy field. 

This same mechanism is responsible for the S3T desta¬ 
bilization of the njc= I perturbation to the jet equilibrium 
at £/£[■= 13.65. However, at e/Sc = 13.65 the jet is un¬ 
stable to both rijc = 0 (with maximum growth rate Or = 
0.324) and to nx= 1 non-zonal perturbations (with maxi¬ 
mum growth rate Or = 0.083 and phase speed Cr = —5.99, 
which is retrograde by 3.18 with respect to the minimum 
velocity of the jet). This equilibrium flow is also hydrody- 
namically stable at all the zonal wavenumbers allowed by 
periodicity (cf. Fig. 5a). This = I unstable eigenfunc¬ 
tion (cf. Fig. 5d) arises from destabilization of the second 
least damped mode, which is the damped external Rossby 
mode indicated in Fig. 5b and shown in Fig. 5c. The en¬ 
ergetics of the instability indicate that the growth of this 
tir = 1 structure arises almost equally from energy trans¬ 
ferred from the mean equilibrium jet to the = I pertur¬ 
bation (aio = 0.160) and energy transferred by the small 
scales (ai> = 0.115) while dissipation accounts for the re¬ 
mainder (7 id = —0.192. 
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Fig. 4. (a) The hydrodynamic stability of at e/Sc — 9. Shown are the maximal modal growth rates, k^Cj, of operator as a function of kx. 
Circles indicate the growth rate at the kx retained in the perturbation dynamics; diamond indicates the growth rate at kx — 1. The equilibrium jet is 
hydrodynamically stable but S3T unstable to = 1 perturbation, (b) The growth rates, kxCj, and phase speeds, c,-, of the least damped eigenvalues 
of A^ for kx — I perturbations. The shaded area indicates the region mm(U^) <Cr< max (U^). The jet is shown in white in both panels (c) and 
(d). The streamfunction of the maximally growing S3T rix — I eigenfunction is shown in (d). This S3T eigenfunction arises from destabilization 
of the least damped mode of A^ with kxCj — —0.15 and c,- = —3.70, indicated with the diamond in (b) and shown in (c). The tix — I S3T instability 
with (7,- = 0.099 and phase speed c,- = —3.81 is supported in this case solely by energy transfer from the mean flow (at the rate CTio = 0.303) 
against the negative energy transfer from the small scale perturbation field (at the rate (7i> = —0.016) and dissipation (at the rate (Tid = —0.188), 
with the growth rate of the S3T instability being Gr — (Jio + + O'id. 
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Fig. 5. (a) The hydrodynamic stability of at e/Sc — 13.65. Shown are the maximal modal growth rates, kxCi, of operator A^ as a function of 
kx. Circles indicate the growth rate at the kx retained in the perturbation dynamics; diamond indicates the growth rate atkx — l. The equilibrium 
jet is hydrodynamically stable but S3T unstable to both fix — 0 and fix — \ perturbations, (b) The growth rates, kxCi, and phase speeds, Cr, of the 
least damped eigenvalues of A^ for kx — \ perturbations. The shaded area indicates the region min((/^) <Cr< max(t/^). The jet is shown in 
white in both panels (c) and (d). The streamfunction of the maximally growing S3T nx—\ eigenfunction is shown in (d). This S3T eigenfunction 
arises from destabilization of the second least damped mode of A^ with = —0.165 and Cr — —6.12, indicated with a diamond in (b) and shown 
in (c). The tix — l S3T instability with G,- — 0.083 and Cr — —5.99 is supported in this case by both energy transfer from the mean flow (at the 
rate aio = 0.160) and energy transfer from the small scale perturbation field (at the rate (7i> =0.115). The dissipation rate is cTid = —0.192. 
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7. Equilibration of the S3T instabilities of the equilib¬ 
rium jet 

We next examine equilibration of the = 1 S3T insta¬ 
bility at e/Cc = 9 and the equilibration of the S3T insta¬ 
bilities at e/Ec = 13.65, which has both = 0 and = 1 
unstable eigenfunctions. 

Consider the energetics of these large scales consisting 
of the kjc = 0 and = 1 Fourier components. Denote 
the kjc = 0 and k^ = 1 components of vorticity of (23a) 
as and and the corresponding vorticity flux diver¬ 
gence of the incoherent components as R^, J?’ and with 

def 

Z*’ = —dyV^ the vorticity of the equilibrium zonal jet. The 
energetics of the equilibration of the S3T instabilities is 
examined by first removing the constant flux to the large 
scales from the small scales that maintains the equilibrium 
flow U‘^. For that reason the vorticity flux divergence asso¬ 
ciated with the deviation of the instantaneous covariance 
from will be considered in the equilibration process. 

Consider first the energetics of the kjc= 1 component of 
the large-scale flow. The first contribution to the energy 
growth of this component is the energy transferred from 
the kx = 0 component of the flow. This occurs at rate; 

<^10 = (Ai„v(t/°)Z',Z1) + (z',Ai„v(C/°)z') , (38) 

with Ainv defined in (34) and t/** the total k^ = Q compo¬ 
nent of the zonal velocity. The second energy source is en¬ 
ergy transferred tokx=\ from the small scales (i.e. those 
with \kx\ > K), which occurs at rate; 

<ri> = (z',/?') + (.Ri,Z‘) , (39) 

with J?' =* —C®) the vorticity flux divergence pro¬ 

duced by covariance C — C*’. Finally, energy is dissipated 
at the rate; 

<riD = (AdZ',Z') + (Z^AdZ*) , (40) 

with Ad defined in (37). 

The energy flowing to the kx = 0 component consist first 
of (§ 01 , the energy transfer rate to this component from the 
kx = I component, which is equal to —§io (being equal 
and opposite to the energy transfer rate to kx = 1 from the 
kx = 0 component), and second of the energy transferred to 
kx = 0hy the small scales, with contribution to the growth 
rate; 

§o> = + (R^\Z^^} , (41) 

with R^ — C*’). Having removed the energy source 

sustaining the equilibrium flow, the energy of Z*’ is dissi¬ 
pated at rate; 

^OD = (Ad(Z'>-Z^),Z'>) + (Z^AojzO-Z^)) . (42) 

The instantaneous rate of change of the energy of the 

and Z* components are then dEo/dt = (§oi + (So> + 


(Sod and dEi/dt = (Sio + (Siy + S\d- By dividing each 
term of dE\/dt with 2(Z*,Z^) we obtain, corresponding 
to (33), (35), (36), the instantaneous growth rates aio, 0'i> 
and aiD and by dividing dEo/dt with 2(Z° — Z®,Z° — Z®) 
the growth rates aoi, Oo> and aoo. As equilibration is ap¬ 
proached the sum of these growth rates approaches zero, 
while the evolution of the growth rates indicates the role 
of each energy transfer rate in producing the equilibration. 

a. Case 1: = 1 instability at e/Ec = 9 

Consider first the equilibration of the «x = 1 instabil¬ 
ity at e/Cc = 9 by first imposing on the jet equilibrium 
the most unstable S3T = 1 eigenfunction at small am¬ 
plitude, in order to initiate its exponential growth phase. 
Evolution of the energy of the Z^ component of the flow 
as a function of time, shown in Fig. 6a, confirms the ac¬ 
curacy of our methods for determining the structure and 
the growth rate of the maximally growing S3T eigenfunc¬ 
tion of the jet equilibrium. The contribution of each of 
the growth rates associated with (38)-(40) to the total nor¬ 
malized energy growth rate of the k,: = \ component of 
the flow, dE\/dt, is shown in Fig. 6b. As discussed ear¬ 
lier, the S3T instability is due to the transfer of energy 
from the zonal flow and the equilibration is seen to be 
achieved by reducing the transfer of the energy from the 
mean flow to the kx = l component by reducing the tilt of 
the non-zonal component of the flow. The Reynolds stress 
contribution remains approximately energetically neutral. 
The flow eventually equilibrates to a nearly zonal config¬ 
uration which is very close to the initial jet, as shown in 
Fig. 7c. The equilibrium state while nearly zonal contains 
an embedded traveling wave (cf. Fig. 7a,b). This wave 
propagates westward with phase speed indistinguishable 
from that of the unstable = 1 S3T eigenfunction, as 
can be seen in the Hovmoller diagram of 'F*, shown in 
Fig. 7d. The PV gradient of the equilibrated jet j3 — d/U^ 
is everywhere positive and the wave propagates in the ret¬ 
rograde part of jet where the PV gradient is close to uni¬ 
form. Also the stmcture of the non-zonal component of 
the equilibrated flow is very close to the structure of the 
most unstable eigenfunction, as seen by comparing Fig. 4d 
with Fig. 7b. This equilibrated state is robustly attract¬ 
ing. When the unstable jet is perturbed with random high 
amplitude perturbations the unstable S3T jet is attracted 
to the same equilibrium. Mixed S3T equilibria of similar 
form have been found as statistical equilibria of the full 
nonlinear equations (Bakas and loannou 2013a, 2014). 

b. Case 2: = 1 instability af e/Cc = 13.65 

The equilibration of the jet at e/Ec = 13.65 involves the 
simultaneous equilibration of two S3T instabilities, of the 
powerful nx = 0 jet instability that grows initially at the 
rate Or = 0.324 and of the weaker n^ = 1 instability that 


Journal of the Atmospheric Sciences 


13 



Fig. 6. (a) Evolution of the disturbance energy, of the deviation of the lai'ge-scale flow from its zonal equilibrium state at e/Cc = 9 with 
equilibrium vorticity Z^. The S3T equilibrium is initially perturbed with the unstable Ux — 1 S3T eigenfunction shown in Fig. 4d. Initially the 
deviation grows at the predicted exponential growth rate of the eigenfunction (dashed) and the equilibration of this instability produces asymptoti¬ 
cally the stationary state shown in Fig. 7a,b comprising a jet with a finite amplitude embedded wave, (b) Evolution of the energetics of the kx — \ 
component of the flow. Shown are the contribution to the instantaneous growth rate of kx = 1 by energy transferred from the mean flow ((7io), from 
the small scales ((Ji>) and that lost to dissipation (cTid)- Also shown is the resulting instantaneous growth rate: (Jir — (Jio + <7i> + (Tid. which 
necessarily vanishes as equilibration is approached. The S3T instability is supported in this case solely from energy transferred to kx — 1 from 
and equilibration is achieved by reducing this transfer. 
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Fig. 7. (a) Mean flow streamfunction, T*, at r = 80, and the velocity field of the kx — 0 and kx — \ components of the equilibrium at e/£c — 9 
resulting from equilibration of the rix — 1 instability. Also shown in white is U^. The equilibrium consists of a jet and a traveling wave that has no 
critical layer in the flow as it travels retrograde with respect to the minimum jet velocity, (b) The wave component of the flow, , and its associated 
velocity field. The wave propagates in the retrograde part of the jet where the potential vorticity gradient, /3 — dyU^ (shown in white) has a small 
and nearly constant positive value, (c) Variation of the zonal flow velocity, U, with y at equilibrium at different x sections. Also shown is 
(dashed line) which is neai'ly identical to the unstable S3T jet U^. (d) Hovmoller diagram of at the location of the minimum of U^, y — 2.4. The 
phase velocity of the equilibrated wave is equal to the phase speed (indicated with the dashed line) of the most unstable S3T eigenfunction shown 
in Fig. 4d. 
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Fig. 8. Evolution of the disturbance energy, 8E,„, associated with the deviation of the large-scale flow Z — Z^, where Z^ is the zonal equilibrium 
vorticity, at e/Sc — 13.65. The S3T equilibrium is initially perturbed with the unstable riy — 0 and rix — I S3T eigenfunctions at small but equal 
amplitude. The tix — 0 eigenfunction grows at Gr — 0.324; the Hx — I at G,- — 0.083 (both indicated with dashed lines). Energy grows at first at the 
rate of the rix — 0 instability, up to f !=s 12, at which time the equilibration of Z^ is established. The equilibration of Z^ is not established until t !=s 60. 
(b) Evolution of the energetics of the Z^. Shown are the contribution to the instantaneous growth rate of Z^ — Z^ from energy transferred: from 
Z^ (c7oi), from the small scales (Ob>)T and that lost to dissipation ((Tod)- Also shown is the actual instantaneous growth rate, (Tor, which vanishes 
at equilibration. The rix — 0 S3T instability is supported by the transfer of energy from the small scales and equilibration is achieved rapidly by 
reducing this transfer, (c) Same as (b) but for Z^. Shown are the transfer rate from Z^ ((Tio), from the small scales (<Ti>) and the energy dissipation 
rate (Tid. The rix — 1 instability is supported by both transfer from Z^ and from small scales and the equilibration is established by reducing the 
transfer from Z^. 



Fig. 9. (a) Mean flow streamfunction, ^,att — 120, and the velocity field of the kx — 0 and^^ = 1 components of the equilibrium at e/Ec — 13.65 
resulting from equilibration of both the fix — 0 and rix — 1 instabilities. Also shown in white is U^. The equilibrium consists of a jet and a traveling 
wave that has no critical layer in the flow as it travels retrograde with respect to all U^. (b) The wave component of the flow, , and its associated 
velocity field. The wave is trapped in the prograde part of the flow where the potential vorticity gradient (shown in white) is large, (c) Variation 
of the zonal flow velocity, U, with y at equilibrium at different x sections. Also shown is (dashed line) which is nearly identical to the unstable 
S3T jet U^. The equilibrated jet is asymmetric, (d) Hovmoller diagram of at the location of the zero of U^,y — 1.6. The phase velocity of the 
equilibrated wave is equal to the phase speed (indicated with the dashed line) of the most unstable S3T eigenfunction shown in Fig. 5d. 
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grows initially at rate Or = 0.083. We impose on the equi¬ 
librium the most unstable S3T = 0 and «x = 1 eigen¬ 
functions at small but equal amplitudes, in order to initiate 
their exponential growth phases. The evolution of the en¬ 
ergy of the Z — Z® component of the flow as a function of 
time (cf. Fig. 8a) shows initial growth at the rate of the 
faster = 0 instability. The equilibration process for the 
n^ = 0 instability is shown in Fig. 8b, and the equilibration 
of the tijc = 1 instability in Fig. 8c. The nx = 0 instability 
is supported by the transfer of energy to the kx = 0 compo¬ 
nent from the small scales (CTo> > 0) as is the equilibrated 
jet. The equilibration of this instability proceeds rapidly 
and is enforced by reduction of the ao>, i.e. the trans¬ 
fer of energy from the small scales. During the equilibra¬ 
tion process there is a pronounced transient enhancement 
of the transfer rate to the mean flow by the eddies. This 
leads to the equilibrated jet shown in Fig. 9a,c which has 
5% greater energy than the original S3T unstable equilib¬ 
rium jet. The equilibrated jet is asymmetric with enhanced 
power at riy = 1. (In this case the unstable «y = 2 jet did 
not merge with the = 1 jet to form a jet with a single jet 
structure.) During the equilibration process aoi is always 
negative, indicating continual transfer of mean jet energy 
supporting the 11^=1 perturbation. The equilibration of 
the % = 1 wave is slower and proceeds in this example, in 
which the jets did not merge, independently of evolution of 
the = 0 instability. The wave is supported by transfer of 
energy from the small scales and from transfer of energy 
from the mean flow. The former remained unaffected dur¬ 
ing the equilibration process and equilibration is achieved 
by reduction of the transfer from the mean flow, aio. The 
PV gradient of the mean flow, j3 — d^U^, shown in Fig. 9b 
is positive almost everywhere and the wave is trapped at 
the prograde part of the jet. As in the case with e/Sc =9, 
the wave propagates at the speed of the S3T eigenfunction 
(cf. Fig. 9d). 

8. Discussion 

a. Correspondence between the S3T dynamics (16) and 

the projected S3T dynamics (23) 

Stability of a two jet state to jet/wave perturbations in 
the projected S3T formulation (23) is shown in Fig. 2. For 
parameters for which the base state becomes unstable to 
non-zonal large-scale perturbations this base state transi¬ 
tions to a new equilibrium in which the jet coexists with a 
coherent wave. The stability calculation, its energetics and 
equilibration process are studied in the framework of the 
projected S3T equations (23), which allows a clear sepa¬ 
ration between the contribution of the coherent jet inter¬ 
action and that of the incoherent eddies to the instability 
and equilibration processes. This stability analysis using 
the projected S3T system produces essentially the same re¬ 
sults as were obtained using the S3T system (16) (compare 
Fig. 2 with Fig. lOa-b and Fig. 3e,f with Fig. lOc-d). The 


equilibrated states produced by these two S3T systems are 
also very similar (cf Fig. 11). 

b. Reflection of ideal S3T dynamics in QL simulations 

The ideal S3T equilibrium jet and jet/wave states that 
we have obtained are imperfectly reflected in single re¬ 
alizations of the flow because fluctuations may obscure 
the underlying S3T equilibrium (cf Farrell and loannou 
(2003, 2014)). The infinite ensemble ideal incorporated in 
the S3T dynamics can be approached in the QL (governed 
by (22)) by introducing in the equation for the coherent 
flow an ensemble mean Reynolds stress obtained from a 
number of independent integrations of the QL eddy equa¬ 
tions with different forcing realizations. 

Consider for example the jet/wave S3T regime at 
e/Cc =9 shown in Fig. 7. The energy of the = 0 
component of the coherent flow is Eq = 1.3 and of the 
kjt= 1 component, which is predominantly aky = l wave, 
is = 0.05. In Fig. 12a-l is shown the approach of the QL 
dynamics to this ideal S3T equilibrium as a function of the 
number of ensemble members, Agns, using as diagnostics 
the structure, indicated by snapshots, of the coherent flow 
and the energy spectrum. Convergence of the energy of the 
QL coherent flow components to that of the S3T as Aens in¬ 
creases is shown in Fig. 12m-n. These ensemble QL simu¬ 
lations were performed by introducing the mean Reynolds 
stress divergence obtained from Aens independent simula¬ 
tions of (22b), all with the same large-scale flow, obtained 
from a single mean QL equation (22a). Convergence to 
the S3T state is closely approached with Aens = 10. In 
simulations with a smaller number of ensemble members 
the ensemble QL supports an irregular weaker ky = 2 jet 
and a stronger kx = l coherent flow, which is concentrated 
at ky = 2 rather than at ky = 1 as predicted in the S3T 
(cf Fig. 12b). As the number of ensemble members in¬ 
creases the jet is more coherently forced and the ideal S3T 
kx = l component, which was previously masked by fluc¬ 
tuations at ky = 2, is revealed. Also note that in these QL 
simulations there are no eddy-eddy interactions and also 
no direct stochastic forcing of the coherent flow compo¬ 
nents and consequently their emergence does not result 
from cascades but from the structural instability mecha¬ 
nisms revealed by S3T. 

Both S3T and ensemble simulations isolate and clearly 
reveal the mechanism by which a portion of the inco¬ 
herent turbulence is systematically organized by large- 
scale waves to enhance the organizing wave. However, 
as in simulation studies revealing this mechanism at work 
in baroclinic turbulence (Cai and Mak 1990; Robinson 
1991), the large-scale wave retains a substantial incoher¬ 
ent component in individual realizations. This is expected 
in the strongly turbulent atmosphere considering that even 
stationary waves at planetary scale which are strongly 
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Fig. 10. (a)-(b): The stability of the jet equilibria in the S3T formulation (16). The corresponding stability properties of the projected S3T 
system are shown in Fig. 2a-b. (c): Contour plot of the streamfunction of the most unstable non-zonal tix — I S3T mean flow eigenfunction of the 
riy — 2 jet equilibrium al e/Sc — 9 with growth rate (7,- = 0.248 and phase speed c,- = — 3.91. The equilibrium jet plotted in solid white. Positive 
(negative) contours are shown with solid (dashed) lines and the zero contour with thick solid line, (b): The energy power spectrum of the mean flow 
eigenfunction. The jet is stable to rix — 0 perturbations but unstable to Hx — I perturbations and the most unstable rix — I eigenfunction is Bloch 
qy — 1 with power atny — 1. 
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Fig. 11. Comparison of the flows resulting from equilibration of the Ux — l S3T instabilities of the projected S3T system (23) (panels (a)-(d)) 
and the S3T system (16) (panels (e)-(h)) at e/Cf = 9. Shown are: The ky energy spectrum of the mean flow {kx — 0) (first column) and of the 
kx — 1 (second column), the kx energy spectrum of both the coherent flow components (|/:;y| < 1, shown with circles) and of the incoherent flow 
components |^j;| > 1 (shown with asterisks) (third row). Snapshots of the mean jet (thick line) and contour plot of the streamfunction of the kx — \ 
wave component are shown in the figures of the fourth column. 
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Fig. 12. (a)-(l): Structure of the mean flow in ensemble QL simulations as the number of ensemble members, Aens^ increases. Shown are 
= 1 (first row, (a)-(d)) and Agns = 10 (second row, (e)-(h)) at e/£c = 9. Also shown is the equilibrated mean flow from the S3T simulation 
(third row, (i)-(l)) at the same forcing amplitude (cf. Fig. 7). Shown are: The ky energy spectrum of the mean flow (kx — 0) (first column) and of the 
kx — 1 (second column), the kx energy spectrum of both the coherent flow components (|^jc| < 1) and of the incoherent flow components |^.v| > 1 
(third row). Snapshots of the mean jet (thick line) and contour plot of the streamfunction of the kx — 1 wave component is shown in the figures 
of the fourth column, (m) The approach of the energy of the zonal flow, Eq , obtained in ensemble QL simulations as the number of ensemble 
members, Agns increases for the case e/Ec — 9. Dashed line marks the S3T prediction, (n): The approach of the wave component energy E\ to 
the S3T predictions (marked by dashed line). Convergence is achieved with Agns = 10. Note that for small Agns fluctuations result in enhanced 
excitation of the kx — 1 component. Parameters as in previous simulations. 
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Fig. 13. QL simulations (panels (a)-(d)), modified NL simulations (panels (e)-(h)) in which only the term (/ — [V • is removed, and 

unmodified ensemble NL simulations (panels (i)-(l)), all with 10 ensemble members at e/gf = 9. Shown are: the ky energy spectrum of the mean 
flow (kx — 0) (first column), of the kx — 1 component (second column), the kx energy spectrum, of both the coherent flow components {\kx\ < 1) 
and of the incoherent flow components kx — 1 (third column). Snapshots of the mean jet (thick line) and contour plot of the streamfunction of the 
wave component = 1 is shown in the figures of the fourth column. Parameters as in the other simulations. Note that the eddy-eddy interactions 
in NL lead in this particular example to an appreciable departure from the QL spectrum of the eddy components, i.e. wavenumbers kx >2 (cf. 
panels (k) and (c)). This is caused because all the eddy components kx > 2 are externally strongly forced (the dissipation has an e-folding of 40 d). 
This figure shows that the eddy-eddy nonlinearity is the source of the difference between the ensemble QL and the ensemble NL simulations. 


forced by topography are revealed clearly only in seasonal 
average ensembles. 


c. Reflection of ideal S3T dynamics in NL simulations 

Consider now the reflection of the S3T jet/wave regime 
in NL and ensemble NL simulations. Ensemble simu¬ 
lations of the NL system (21) were performed by intro¬ 
ducing in the mean equations (21a) the ensemble aver¬ 
age of Pk (u' • Vi^') and (U • V -I- u' • VZ) obtained 
from A^ens independent simulations of the perturbation 
NL equations (21b) all with the same large-scale flow. 
The corresponding results of the ensemble QL simulation 
(cf. Fig. 13a-d) differ from those of the ensemble NL sim¬ 
ulation. The nonlinear term (/ — Pk) (u' • Vi^') is respon¬ 
sible for the difference between the NL and QL ensemble 
simulations, as shown in Fig. 13e-h. In this figure an en¬ 
semble integration of the NL equations with this term ab¬ 
sent is shown to produce results that are very close to the 
QL results. 

When all waves with \kf\>2 forced equally, as in 
the S3T examples discussed above, the eddy-eddy inter¬ 
actions are strong in the corresponding NL resulting in a 


substantial modification of the spectrum of the eddy mo¬ 
tions which is not reflected in the S3T. To obtain corre¬ 
spondence an effective stochastic forcing which param¬ 
eterizes the absent eddy-eddy interactions is required in 
S3T (Constantinou et al. 2014). 

Alternatively, when the term (7 —Pa") [V • (u'^')] is sup¬ 
pressed by choosing low forcing excitation, which results 
in weak modification of the spectrum of the incoherent 
component, agreement between NL and QL simulations 
is obtained. This is demonstrated in Fig. 14, in which 
we show results obtained with an approximate small-scale 


if Kf -6Kf<k< Kf -b 6Kf , 


if \k-Kf \ > 5KfOi- \k^\ < 1 , 
(43) 

with Kf = 10, 5Kf = 1 and r = 0.01, v = 0 , as in Bakas 
and loannou (2014). With these parameters the S3T zonal 
jet equilibrium is stable to jet perturbations and unstable 
to = 1 wave perturbations and the resulting equilibrium 
state in NL has a wave k^ = 1 component in agreement 
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Fig. 14. Comparison among QL and NL simulations with ensemble members A/ens = 1 (panels (a)-(h)), A^ens = 80 (panels (i)-(p)) and S3T 
simulation (panels (q)-(t)). The S3T state in this example predicts a jet component and a = 1 component and this is reflected in both QL and 
NL ensemble simulations as revealed by the ky energy spectrum of the respective mean flow {kx — 0) (first column), the kx — l component (second 
column), and the kx energy spectrum (third column). Snapshots of the mean flow (thick line) and contour plot of the streamfunction of the wave 
component kx — l are shown in panels (d), (h), (1), (p) and (t) (fourth column). Simulations with isotropic forcing at total wavenumber Kf — 10 
(cf. (43)) with e — 4.2 x 10“^, linear damping coefficient r — 0.01 and v — 0. 
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Fig. 15. Comparison among QL and NL simulations with ensemble members A/ens = 1 (panels (a)-(h)), A^ens = 80 (panels (i)-(p)) and S3T 
simulation (panels (q)-(t)). The S3T state in this example predicts only a jet component and no kx — 1 component and this is reflected in both QL 
and NL ensemble simulations as revealed by the ky energy spectrum of the respective mean flow {kx — 0) (first column), the kx — I component 
(second column), and the kx energy spectrum (third column). Snapshots of the mean flow (thick line) and contour plot of the streamfunction of 
the wave component kx — \ are shown in panels (d), (h), (1), (p) and (t) (fourth column). Simulations with the anisotropic forcing spectrum (25) at 
e — 0.830 with zonal wavenumbers 8 < \kx\ < 14 being forced. Mean linear damping coefficient is rm = 0.1 and linear damping coefficient of the 
incoherent flow is r = 1.5 and v = 0. 
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with S3T predictions. Also the energetics of the mecha¬ 
nism of destabilization of the external = 1 Rossby wave 
is partitioned between coherent and incoherent sources 
consistent with the mechanism described in the previous 
sections. 

It could be maintained that because isotropic ring forc¬ 
ing suppresses eddy-eddy interactions, the agreement be¬ 
tween S3T and NL should be expected (cf. Bakas et al. 
(2015), Appendix C). This property follows from the fact 
that a barotropic fluid excited in an infinite channel with an 
isotropic ring forcing with spectrum 2(k) 5(|k| —Kf) 

results in a nonlinear solution which by itself could never 
give rise to a jet. The emergence of jets under this forcing 
can only result from imposition of a separate perturbation 
such as the jet perturbation that results in the S3T jet in¬ 
stability. As an example closer to physical reality consider 
forcing of the form (25) which excites the zonal wavenum¬ 
bers 8 < \kx\ < 14, modeling baroclinic energy injection, 
and with linear damping for the incoherent scales r = 1.5 
and Tm = 0.1 for the coherent scales with corresponding 
e-folding times of 4d and 60d. For these parameters S3T 
theory predicts that the 5 jet equilibrium at e = 0.830 is 
stable to both jet and = 1 wave perturbations and con¬ 
sequently S3T theory predicts that QL and NL simulations 
should show suppressed energy in the = 1 coherent 
wave component. Good agreement between QL and NL 
is found in the channel as shown in Fig. 15. It is inter¬ 
esting to note that while the jet in the NL simulation has 
the structure predicted by the S3T its amplitude is reduced 
consistent with a component of the eddy-eddy interactions 
acting as diffusion on the mean jet. 

9. Conclusions 

Large-scale coherent structures such as jets, meander¬ 
ing jets and waves embedded in jets are characteristic fea¬ 
tures of turbulence in planetary atmospheres. While con¬ 
servation of energy and enstrophy in inviscid 2D turbu¬ 
lence predicts spectral evolution leading to concentration 
of energy at large scales, these considerations cannot pre¬ 
dict the phase of the spectral components and therefore 
can not address the central question of the organization 
of the energy into specific structures such as jets and the 
coherent component of planetary scale waves. In order 
to study structure formation additional aspects of the dy¬ 
namics beyond conservation principles must be incorpo¬ 
rated in the analysis. For this purpose SSD models have 
been developed and used to study the formation of coher¬ 
ent stmcture in planetary scale turbulence. In this work 
an SSD model was formulated for the purpose of study¬ 
ing the regime of coexisting jets and waves. In this model 
a separation in zonal Fourier modes is made by projection 
in order to separate a coherent structure equation, in which 
only the gravest zonal harmonics are retained, from the re¬ 
maining spectrum, which is assumed to be incoherent and 


gives rise to the ensemble mean second order statistics as¬ 
sociated with the incoherent turbulence. This second order 
SSD model is closed by a stochastic forcing parameteriza¬ 
tion that accounts for both the neglected nonlinear dynam¬ 
ics of the small scales as well as the forcing maintaining 
the turbulence. The equation for the large scales retains 
its nonlinearity and its interaction through Reynolds stress 
with the perturbations. 

In this model jets form as instabilities and equilibrate 
nonlinearly at finite amplitude. A stable mode of the 
Rossby wave spectrum associated with these jets is desta¬ 
bilized for sufficiently strong forcing by interaction with 
perturbation Reynolds stresses. This destabilization is 
found to have in some cases the remarkable property of re¬ 
sulting from destabilization of the retrograde Rossby wave 
to mean jet interaction by structural modification of this 
damped mode arising from its interaction with the inco¬ 
herent turbulence thereby transforming it into an unstable 
mode of the mean jet. In other cases comparable con¬ 
tributions are found from direct forcing by the Reynolds 
stresses, as in S3T instability with projections at = 0, 
and induced jet/wave interaction, as in traditional hydro- 
dynamic instability. This synergistic interaction provides a 
powerful new mechanism for maintaining planetary waves 
that will be the subject of further investigation. 
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